Escape from the potential well: competition between long jumps and long waiting 
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Within a concept of the fractional diffusion equation and subordination, the paper examines the 
influence of a competition between long waiting times and long jumps on the escape from the poten- 
tial well. Applying analytical arguments and numerical methods, we demonstrate that the presence 
of long waiting times distributed according to a power-law distribution with a diverging mean leads 
to very general asymptotic properties of the survival probability. The observed survival probability 
asymptotically decays like a power-law whose form is not affected by the value of the exponent char- 
acterizing the power-law jump length distribution, ft is demonstrated that this behavior is typical 
of and generic for systems exhibiting long waiting times. We also show that the survival probability 
has a universal character not only asymptotically but also at small times. Finally, it is indicated 
which properties of the first passage time density are sensitive to the exact value of the exponent 
characterizing the jump length distribution. 
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I. INTRODUCTION 

Kramers' seminal paper [l[ studies the problem of the 
escape of a particle over the potential barrier. An escape 
of the test particle from the metastable state takes place 
due to the presence of noise. Otherwise, in the static 
potential, any particle (described by the Newtonian me- 
chanics) would be sliding to the nearest minimum of the 
potential only. Random interactions of the test particle 
with other particles of the medium make the escape from 
the metastable state possible. 

In the simplest situations, it is assumed that the noise 
acting on a test particle is white and Gaussian. In far 
from equilibrium situations, the stochastic process de- 
scribing interactions of the test particle with the envi- 
ronment can be still of the white type but its values can 
be distributed according to some heavy tailed distribu- 
tion. An a-stable noise, leading to Levy nights [HI], 
can be considered 0-0] as an exemplary driving process. 
Levy stable noise is especially suitable for description of 
the far from equilibrium situations because of its heavy- 
tail characteristics, infinite divisibility and the presence 
of the generalized central limit theorem 0] ■ The presence 
of cv-stable Levy type noises affects the escape process. 
However, the survival probability, as in the case of the 
white Gaussian noise, still has the exponential asymp- 
totic form P-Q. Levy flights which have been observed 
in various real-life situations when a test particle is able 
to perform unusually large jumps pi. ficl ITll] . Levy flights 
have been documented to describe motion of fluorescent 
probes in living polymers, tracer particles in rotating 
flows [l2j and cooled atoms in laser fields. They also 
serve as a paradigm of efficient searching strategies on 



folding polymer chains 13| and environmental problems 

El- 

The deviations from the exponential asymptotic form 
of the survival probability can be introduced by the trap- 
ping events distributed according to a power- law dis- 
tribution with the diverging mean [TJ [l5j], see Sec. [TTJ 
Transport in porous, fractal-like media or relaxation ki- 
netics in inhomogeneous materials lead to long trapping 
events which can be distributed according to a power- 
law. Consequently, such processes are usually ultraslow, 
i.e. subdiffusive [Tl| . The most intriguing situation takes 
place, however, when both effects the presence of long 
waiting times for the next step and l ong jum ps - are com- 
bined in the same scenario, see fTTL fl6P4]~9| . because long 
jumps can hide temporal dynamics induced by long wait- 
ing times. This competition can be of special relevance 
in single particle tracking experiments [20j. For example, 
a particle moving along a polymer can be trapped in a 
series of local minima of the potential while long jumps 
can take pla ce in conformational space of the fast folding 
polymer |13l [2l| where the distance measured along the 
polymer can be significantly longer than the Euclidean 
distance. In such a case instead of travelling a long dis- 
tance along the polymer, due to chain looping and volume 
exchange [13| |. it is possible to jump directly to the other 
segment of the polymer. 

In the region of non-Markovian Levy flights [ll|, i.e. 
continuous time random walks with power-law waiting 
time p(At) oc t~( v+r) (0 < v < 1) and jump length 
p(Ax) ex |a;|~( Q+1 ) (0 < a < 2) distributions , the evolu- 
tion of the probability density of finding a random walker 
at the time t in the vicinity of x is asymptotically de- 
scribed by the fractional diffusion equation [n], l22j 
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where V(x) represents external, static potential. In 
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Eq. ((T|), d a /d\x\ a stands for the Riesz fractional (space) 
derivative which is defined via the Fourier transform 
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Here, we focus on the investigation of the influence of a 
competition between long waiting times and long jumps 
on the Kramers problem [l|, l24j - |2(| . As a jump length 
distribution we chose symmetric a-stable densities which 
are responsible for the occurrence of the fractional space 
derivative in Eq. ((TJ) [27], [H[ . Symmetric a-stable den- 
sities asymptotically behave like a power- law, see @-[|j|. 
Nevertheless, what is even more important is the fact 
that for jump lengths p(Ax) distributed according to 
symmetric a-stable densities the CTRW scenario leads 
exactly to Eq. (UJ). The parameter a (0 < a ^ 2) is called 
the stability index, see [2J . In the limit of a = 2 any a- 
stable density is equivalent to the Gaussian distribution. 
Therefore, in the limit of a = 2 the fractional (space) 
derivative in Eq. ((TJ) is replaced by the partial deriva- 
tive of the second order. The parameter v (0 < v < 1), 
see Eq. (UJ), is responsible for the trapping events and is 
traditionally called the subdiffusion parameter [ll|. In 
the limit of v = 1 the process described by Eq. (JTJ) be- 
comes a Markovian diffusion process and the fractional 
(time) derivative disappears. Breaking of Markovianity 
is introduced by the long waiting times which are de- 
scribed by the fractional (time) derivative [TJ [HI- Fi- 
nally, in the limit of a > 2 and v > 1 the mean waiting 
time and variance of the jump length becomes finite and 
the fractional Fokkcr-Planck equation is replaced by the 
standard Fokker-Planck equation (29l - l32j . 
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FIG. 1. The Kramers problem is studied for the generic 
double-well, fourth-order potential V(x) = — fa; 2 + \x 4 . Ini- 
tially, a test particle is located in the left potential well. 



II. METHODS AND RESULTS 

We extend earlier studies [!, 0, HH by investigating 
the influence of the competition between long waiting 



times and long jumps on the Kramers problem [l|, [23 - [2r| . 
Initially, at the time t = a test particle is located in the 
left potential well of the generic, fourth-order, double- 
well potential 



V(x) 



(2) 



The minima of the potential, see Eq. ([2]) and Fig. [TJ are 
located at ±x w = ±^/a/6. The height of the potential 
barrier is AV = V(0) - V(x w ) = a 4 /4b. We study and 
compare two different setups: (i) escape from the poten- 
tial well to the barrier top (x^) and (ii) escape from the 
potential well to the (other) potential well (x w ) Q, i.e. 
we assume that the border between the states is located 
either at the barrier top or at the other potential well, 
see Fig. [TJ 

In the case of (Markovian) Levy flights (y = 1 with 
a < 2) Eq. (UJ) reduces to 



dp(x, t) 
dt 



d d a 



p(x,t). 



(3) 



The stochastic representation of solutions of Eq. §3§ is 
provided by the following Langcvin equation 



dX(t) = -V'(X(t))dt + dL afi {t), 



(4) 



driven by the symmetric a-stable motion L a _o(t) whose 
increments, AL at o(At) = L a $(t + At) — L a ${t), are 
independent and distributed according to the symmet- 
ric a-stable density with the corresponding characteris- 
tic function 4>{k) = (cxp(ifcAL Q ,.o)) = exp (— At<r a \k\ a ) 
0-13l- Here the parameter a € (0, 2] denotes the stability 
index, describing the asymptotic long-tail power-law for 
the ^-distribution, which for a < 2 is of the |x| _ ' 1+0! ' ) 
type. The parameter a (a € (0,oo)) is the scale param- 
eter which characterizes the overall distribution width. 

For the exponent characterizing the waiting time dis- 
tribution v < 1, the stochastic representation of solu- 
tions of Eq. (UJ) is provided by the subordination method 
HH |35{ . i.e. the process inspected X(t) is obtained as 
a function X(t) = X(S t ) by using a different (stochas- 
tic) clock St- In this framework St denotes the zz-stable 
subordinator, S t = inf {s : U(s) > t}, where U(s) stands 
for a strictly increasing z^-stablc process whose distribu- 
tion has the Laplace transform (e~ ku ^) = e~ sk . The 
parent process X(s) is composed of independent incre- 
ments of the symmetric a-stable motion described in an 
operational time s 



dX{s) = -V'(X(s))ds + dL a . {s). 



(5) 



The subordination provides a proper stochastic realiza- 
tion of the random process described by the fractional 
Fokker-Planck equation (UJ) [13, HH • The stochastic rep- 
resentation of solutions of Eq. (UJ), see [13, HU, is valid 
not only asymptotically (large x and long t) but for all 
times t and positions x [34l . |35| . The stochastic repre- 
sentation of solutions of Eq. (UJ) based on the subordi- 
nation method does not use explicitly the waiting time 
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distribution. The information about waiting time distri- 
bution is implicitly incorporated in the ^-stable subordi- 
nator St which provides a link between the physical time 
t and the operational time s [HI, (3||. However, there 
are other, asymptotic, stochastic representation of solu- 
tions of Eq. ([TJ which make explicit use of the waiting 
time distribution [3IH Wf\ , when indeed the waiting time 
distribution has the power-law asymptotics. 

The first passage time density is related to the proba- 
bility density p(x,t), see Eq. (P) 



= -| £j(x,t)dx = -±S^(t), 



(6) 



where S w ~ x (t) = J^ oo p(x,t)dx is the survival probabil- 
ity, i.e. the probability that at the time t a test particle, 
which started its motion at the left potential well, has not 
been yet absorbed at the boundary located at the barrier 
top Xh or at the potential well x w . In more general cases, 
when trajectories become discontinuous, boundary con- 
ditions are non-local, see below. The mean first passage 
time can be calculated as 



rriVJ — X 

MFPT 



tr~*{t)dt 



S w -*(t)dt, (7) 



where x stands for the final point, i.e. the barrier top 
(xb) or the potential well (x w ). In order to guarantee 
the finite value of the mean first passage time, the sur- 
vival probability (as well as the first passage time den- 
sity) asymptotically has to tend to zero fast enough. In 
particular, for v = 1 and a = 2, the standard Kramers 
problem [l| is recovered, the survival probability has the 
exponential asymptotics and the mean first passage time 
can be calculated exactly [13]. The survival probability 
attains the exponential asymptotics 0-01 also for Marko- 
vian Levy flights when the mean first passage time is fi- 
nite [HQ- However, due to memory effects for v < 1 the 
asymptotic behavior of the survival probability is of the 
power-law type, see below. In such a case, in order to 
assure a finite value of the mean first passage time, the 
exponent characterizing the decay of the survival prob- 
ability needs to be larger than one, i.e. S w ~*(t) cx 
with 7 > 1. By analogy, the decay of the first passage 
time density needs to be characterized by the exponent 
larger than two. Otherwise, the mean first passage time 
would be infinite due to divergence at infinity. 

Properties of the first passage time distribution can be 
determined from Eqs. ([1]) an d (l6l) . The solution p(x,t) 
of Eq. ([T]) can be written [ll], I30L [32j as a sum of eigen- 
functions 



POM) = ^2c l p i {x,t) =^2c t T t (t)if l (x) 



(8) 



where T^t) fulfills 

dTj{t) 
dt 



(9) 



while ifi(x) fulfills 



dx y ' d\x\ a 



ipi(x) = -\i, u tpi(x). (10) 



The solution of Eq. © is given by the Mittag-Lcfflcr 
function 12311 



T i {t)=E v {-\ itV t v ) = Y j 



r(i + vj) 



(ii) 



Eigenvalues A^ are defined by boundary conditions to 
Eq. (|10p and coefficients c» are determined by the initial 
condition. More precisely, at the barrier top (xb) or at 
the potential well (x w ) there are absorbing boundaries. 
Consequently, for a = 2 one has p(x — Xb,t) = or 
p(x = x w , t) — while for a < 2 (due to discontinuity of 
trajectories) onehasp(x X\,,t) — 0orp(x x w ,t) = 0, 
see [H. 

The main features of the first passage time distribution 
can be deducted without the knowledge of the eigenval- 
ues \i tV because of the properties of the Mittag-Lcfflcr 
function. The Mittag-Lcfflcr function has a stretched ex- 
ponential form for small values of the argument 



E v [-t v ) oc exp 



(12) 



while asymptotically for large arguments it behaves like 



E v (-t u ) (x r 



(13) 



llj . For v < 1, the asymptotic behavior of the 



Mittag-Lcfflcr function is solely determined by the value 
of the exponent v. In the limit of v = 1 the Mittag-Lcfflcr 
function is equivalent to the exponential function. 

The survival probability, S(t), and the first passage 
time density, /(i), can be calculated by use of Eqs. (|6]) 
and ((8]). The integration of the spatial eigenfunctions 
fi(x) can be considered as a modification of constants 
Ci, see Eq. ([5]). Therefore, S(t) is a sum (with unknown 
coefficients) of the temporal eigenfunctions Ti(t). For 
v = 1, the process studied is Markovian and the sur- 
vival probability has the asymptotic exponential behav- 
ior which is determined by the (unknown) smallest eigen- 



value A min = min^Ai} 

S„=i(t) oc exp(-A mi „f). 



(14) 



For v < 1, the process studied is non-Markovian and 
the survival probability has the asymptotic power-law 
behavior (sec Eq. (fl"3"|) ) determined by the exponent v 



S(t) oc r 



(15) 



Figures present numerical results obtained using 
the subordinated Langevin equation, see Eq. ([S]) and 
pjil [35| . Numerical simulations were performed with the 
same set of parameters as in 0], i.e. a = y/2, a = 128 
and b = 512 (see Eq. @ and Fig. [T])). Such a choice of 
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parameters leads to x w = 1/2 and AV = 8, i.e. the po- 
tential wells are located at ±1/2 while the barrier height 
is 8. Numerical simulations have been performed with 
the time step of integration At = 10 -4 and averaged 
over N = 3 x 10 5 realizations. The range of parameters 
has been restricted to 0.5 ^ a ^ 2 and 0.5 ^ v ^ 1, 
because a small value of the stability index a might lead 
to numerical inaccuracy while small values of the wait- 
ing time exponent result in a significant increase of the 
simulation time. In the limit of v = 1, numerical meth- 
ods applied perfectly recovered all the results of earlier 
studies @. Therefore, the Monte Carlo inspection of the 
Markovian Kramers problem served as a benchmark for 
implemented numerical methods. 




t t 



FIG. 2. The rescaled survival probability S(t) for the Kramers 
problem in the case of non-Markovian Levy flights. The ab- 
sorbing boundary is located at the barrier top (t) - 
left panel) or at the right potential well x w (S w ~ w (t) - right 
panel). The value of the waiting time exponent is v = 0.95. 
Various symbols correspond to various values of the stability 
index a. The solid lines present the survival probability for 
the escape from the finite interval of the length L = 1, see 




In order to estimate the survival probability it is 
necessary to collect the sample of first passage times. 
Required first passage times are obtained by numeri- 
cal simulation of the subordinated Langevin equation 
SHESIllim- Th e Langevin equation © is simu- 
lated until the first escape from the domain of motion, i.e. 
as long as x(t) < x^ or x(t) < x w . For a < 2, trajecto- 
ries of the process x(t) are discontinuous. Consequently, 
boundary conditions become non-local [H| because the 
barrier location is not hit by a majority of trajectories. 
Therefore, in order to assure correct evaluation of the first 
passage time it is necessary to impose non-local bound- 
ary conditions. Within simulations it is assumed that the 
entire semi line x ^ Xh or x ^ x w is absorbing. The dis- 
continuity of trajectories is responsible for failure of the 
method of images [ill, E3 an d the presence of leapovers 
[i^ . Due to leapovers, the particle escaping from the ini- 
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FIG. 3. The same as in Fig. [2] for the waiting time exponent 
v — 0.7. The asymptotic power-law behavior, characterized 
by exponent — v, is well visible. 



tial potential well can jump to a distant point located far 
beyond the position of the absorbing boundary. With the 
non-negligible probability, the jump can be long enough 
to overpass simultaneously both locations of absorbing 
boundaries Xb and x w . 

Figures [2] and [3] present survival probabilities for a sys- 
tem in which there is a competition between long wait- 
ing times and long jumps. The survival probabilities de- 
picted in Figs. [5] and [3] have been rescaled in such a way 
that medians of first passage times densities are located 
at the same point, see [3!|. The insets show the behavior 
of the survival probability at small values of first passage 
times. Furthermore, the solid lines present the depen- 
dence of the survival probability for the escape from the 
finite interval of width L = 1 , see [39j , whose asymptotics 
is of the same type as the asymptotic dependence of the 
survival probability in the generalized Kramers problem. 
In the case of both the anomalous diffusion on finite inter- 
vals and the Kramers problem the asymptotic power-law 
behavior is determined by the waiting time exponent v. 
Figs. [2] and [3] confirm that the asymptotic dependence of 
survival probabilities is the same regardless of the exact 
value of the stability index a, see Eq. (fT5|) . The survival 
probability from the finite intervals and survival proba- 
bilities for the Kramers problem differ at small values of 
first passage times, see insets in Figs. [2] and [3] In the limit 
of the waiting time exponent v = 1, the studied process 
becomes a Markovian Levy flight and the survival proba- 
bility attains its exponential asymptotic dependence, see 
Eq. CL1 and jjjjji. 

The asymptotic power-law decay of the survival prob- 
ability S(t), sec Eq. (15]) and Figs. [2] and |3] is the di- 
rect consequence of the properties of the Mittag-Lcfflcr 
function which describes decay of single modes of the 
fractional Fokker-Planck equation (JTJ). The power-law 
decay observed in non-Markovian situations (v < 1) can 
be contrasted with Markovian Levy flights (y = 1) asso- 
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ciated with the exponential asymptotics of the survival 
probability. Moreover, non-Markovian properties of the 
system studied are not only demonstrated by a slow de- 
cay of the survival probability but also by a long mem- 
ory of initial conditions [HI, [3^, H|| . Both long memory 
signatures mentioned originate in temporal properties of 
the fractional Fokker-Planck equation. In other words, 
non-Markovian Levy flights (anomalously long jumps co- 
existing with long waiting times) do not diminish the 
memory effects of the process studied. The presence of 
long jumps can be observed in the tail asymptotics of 
p(x,t), see [i~8l l45j. In such a situation the waiting time 
exponent v determines the ratio of convergence to the 
asymptotic dependence or to stationary states [IB)]. Fi- 
nally, the slowly decaying memory (controlled by the ex- 
ponent v) is responsible for a slow asymptotic decay of 
the survival probability (characterized by the exponent 
v < 1) and the divergence of the mean first passa ge t ime 
(as in the case of escape from finite intervals [39|, |46| ) . 

Figures 2] and [5] present first passage time densities for 
a small value of first passage times for values of the wait- 
ing time exponent v = 0.95 (Fig. [4]) and v = 0.7 (Fig. [5]). 
First passage time densities are rescaled in the analogous 
manner as survival probabilities. The rescaled densities, 
for both the potential well to the barrier top (f w ~ h (t)) 
and the potential well to the potential well (/ w ~ w (i)) 
setups, are in the perfect agreement, i.e. curves corre- 
sponding to various values of the stability index a are 
not distinguishable. For small and moderate values of 
the waiting time exponent, the short time behavior of 
the first passage time densities agrees with the short 
time asymptotics derived from the short time form of the 
Mittag-Lcffler function (solid lines), see Eq. (fT2)) . How- 
ever, the level of agreement is better for the potential well 
to the barrier top setup than for the potential well to the 
potential well configuration. Furthermore, with the de- 
crease of the waiting time exponent the level of agreement 
increases, see Figs. [4] and [5j because a decreasing value 
of the waiting time exponent increases the probability of 
long waiting times to occur. 

In a general case, the asymptotic of the survival proba- 
bility is determined by the value of the waiting time expo- 
nent v and properties of the Mittag-Lcffler functions. For 
large values of the first passage time t, survival probabil- 
ities behave like t~ v {v < 1), see also [15|]. Consequently, 
analogously as for the escape from the finite intervals 
[39l l4q , the mean first passage time from the potential 
well diverges. Due to a slow power-law decay of the first 
passage time density it cannot be characterized by cither 
the mean value or the variance. Therefore, in order to 
further elucidate properties of the first passage time den- 
sity we investigate properties of the median and width of 
the first passage time distributions. The median is the 
quantile (70.5 of the first passage time density, i.e. it is 
defined by J^ ' 5 f(t)dt = 0.5 or equivalent S^go.s) = 0.5. 
To characterize the width of the distribution we use the 
interquantile distance go. 9 — Q0.1, he. the width of the 
interval containing 80% of first passage times. 
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FIG. 4. The rescaled first passage time density f(t) for 
the Kramers problem in the case of non-Markovian Levy 
flights. The absorbing boundary is located at the barrier top 
x b (/ w ~ b (£) ~~ left panel) or at the right potential well x w 
(/ w ~ w (t) - right panel). The value of the waiting time ex- 
ponent is v — 0.95. Various symbols correspond to various 
values of the stability index a. The solid lines present short 
time behavior of the derivative of the Mittag-Leffler function. 
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FIG. 5. The same as in Fig. [4] for the waiting time exponent 
v = 0.7. 



Figure [6] presents the location of the median as a func- 
tion of the stability index a for the potential well to the 
barrier top (q^ J b - left panel) and the potential well to 
the potential well (9^5 w ~ right panel) setups. Fig. [7] 
presents the width of the first passage time density as a 
function of the stability index a for the potential well to 
the barrier top (^9 b — 1oi h ~ panel) and the poten- 



tial well to the potential well (q™ 



w— w 

lo.l 



9 -H0.1 ri g nt panel) 
setups. Analogously, Figs. [8] and [9] present the location 
of the median and the distribution width as a function 
of the waiting time exponent v. 

Both the location of the median (90.5) and the distri- 
bution width (<7o.9 — <7o.i) are non-monotonic functions of 



6 



2 | 1 1 1 2 




0.4 1 1 1 1 0.4 1 1 1 1 

0.5 1 1.5 2 0.5 1 1.5 2 

a a 



FIG. 6. Locations of the medians (quantile go. 5) of the first 
passage time densities as the function of the stability index a 
for the potential well to the barrier top {q^^ b - left panel) and 
the potential well to the potential well (5*5" - right panel) 
setups. Various symbols correspond to various values of the 
waiting time exponent v. Lines are drawn to guide the eye 
only. 
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a a 

FIG. 7. The width of the first passage time density defined 
as the interquantile distance (go. 9 — qo.i) as the function of 
the stability index a for the potential well to the barrier top 
(l™9 b ~ 1oi h ~ l e ft panel) and the potential well to the poten- 
tial well (q$g w — q^i™ - right panel) setups. Various symbols 
correspond to various values of the waiting time exponent v. 
Lines are drawn to guide the eye only. 



the stability index a with the minima located at a w 1.5, 
see Figs. [6] and [7] This effect might be considered as the 
analogy to the fact that in the presence of (Markovian) 
Levy flights (v = 1) the mean first passage time is a non- 
monotonic function of the stability index a @ . Due to 
long waiting times the studied process is characterized 
by the infinite mean first passage time, however, the me- 
dian (go. 5) as a function of the stability index a behaves 
qualitatively in the same way as the mean first passage 



time for v = 1 @. For the potential well to the barrier 
top setup, the minima of the median and the distribution 
width are located at slightly larger values of the stability 
index a than for the potential well to the (other) poten- 
tial well setup, see Figs. [6] and [7] Finally, in the limit of a 
small value of the stability index a (with the fixed value 
of the exponent v) both setups become indistinguishable. 

The dependence of the location of the median ((70.5) 
and the distribution width ((70.9 ~ <7o.i) as the function of 
the waiting time exponent v is more regular, see Figs. [5] 
and [9] More precisely, locations of the median exhibit 
some weakly non-monotonic behavior, see bottom curves 
in Fig. [5J while the distribution width decreases with the 
increase of the exponent v, see Fig. [5] The dependence 
of the width of the first passage time distribution as a 
function of the waiting time exponent is typical, i.e. the 
width of the distribution decreases with the increase of 
the exponent v. This monotonic dependence can be ex- 
plained by the fact that with increasing i/, long waiting 
times for the next step become less probable, which in 
turn leads to a decrease in the fraction of extreme events 
and decrease of the distribution width. 

Despite the fact that for v < 1 the mean first passage 
times (?Mpp T and ?Mpp T ) diverge, it is still possible to 
find an optimal value of the stability index a. The op- 
timal value of the stability index a leads to the minimal 
value of the median (see Fig. [SJ or the minimal width 
of the first passage time distributions (see Figs. [7]) • This 
non-monotonic dependence of median locations or distri- 
butions' widths can be explained by the non-equilibrium 
character of Levy noises [2l], 53 ■ On the one hand, fol- 
lowing the line of reasoning sketched in [2l[ , fluctuations 
in temperature in the Langevin equation driven by the 
white Gaussian noise (£(t)), i.e. ^ = y / 2T(t)£ > (t), re- 
sult in the occurrence of Levy noises. In other words, if 
fluctuations in temperature T(t) are ruled by a one sided 
Levy process £ a /2,i(i)j it is possible to derive the frac- 
tional Fokker-Planck equation (JXJ) with the spatial frac- 
tional derivative using the Langevin equation with the 
fluctuating temperature [2l|. On the other hand, fluc- 
tuations in temperature alter the height of the potential 
barrier, because within equilibrium statistical mechan- 
ics the height of the potential barrier can be related to 
temperature. Therefore, the system studied resembles 
a classical model of the resonant activation [48[ and the 
mean first passage time (for v = 1), as well as the median 
and the width of the first passage time distribution are 
non-monotonic function of the stability index a. Never- 
theless, this analogy is very loose. On the contrary to the 
Gaussian white noise, due to the non-equilibrium char- 
acter of white Levy noises (with a < 2) the distribution 
width is not related to temperature. 



III. CONCLUSIONS 

Using the recently developed scheme of the subordina- 
tion (34l HH, we have studied the Kramers problem in 
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FIG. 8. The same as in Fig. [6] as the function of the waiting 
time exponent v. Various symbols correspond to various val- 
ues of the stability index a. Lines are drawn to guide the eye 
only. 
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FIG. 9. The same as in Fig. [7|as the function of the waiting 
time exponent v. Various symbols correspond to various val- 
ues of the stability index a. Lines are drawn to guide the eye 
only. 



the case of non-Markovian Levy flights, i.e. long waiting 
times coexisting with long jumps distributed according to 
a-stable densities. The model considered exhibits a be- 
havior similar to escape from the positive semi-line [49i — 
[5f| and finite intervals [39j when the asymptotic behavior 
of the first passage time distribution is not sensitive to 



the details of the jump length distribution but is deter- 
mined by the waiting time exponent v. 

Using analytical and numerical arguments, we have 
demonstrated that the asymptotic power-law behavior 
of the survival probability is a general feature of the 
non-Markovian Levy flights which are characterized by 
anomalously long waiting times distributed according 
to a probability density with a power-law asymptotics 
and jump lengths distributed according to the a-stable 
density. Analytical arguments are provided by proper- 
ties of the Riemann-Liouville fractional time derivative 
which leads to the occurrence of the Mittag-Lcffier de- 
cay pattern, while the numerical arguments are based 
on the Monte Carlo simulations of trajectories of the 
stochastic process otherwise described by the fractional 
Fokkcr-Planck equation. The numerical results con- 
structed by the subordination method (which provides 
a stochastic representation of solutions of the fractional 
Fokker-Planck equation) corroborated analytical predic- 
tions. Furthermore, by a proper rescaling of first pas- 
sage time densities, we have demonstrated their universal 
character not only for large values of first passage times 
but also for short times. 

For the fractional Kramers problem, the mean first pas- 
sage time characterizing the escape of the particle over 
the potential barrier is infinite. The survival probability 
is asymptotically characterized by the exponent which is 
equal to the value of the exponent characterizing asymp- 
totic power-law decay of the waiting time distributions. 
The power-law decay of the survival probability is not 
fast enough to ensure a finite value of the mean first pas- 
sage time. The presence of long jumps (coexisting with 
long waiting times) does not change the asymptotic be- 
havior of the survival probability. Consequently, despite 
the presence of long jumps, the mean first passage time 
for the fractional Kramers problem remains infinite. 

For the value of the waiting time exponent v < 1, 
the mean first passage time cannot be used to charac- 
terize the fractional Kramers problem. The escape of 
the particle from the potential well can be quantified by 
properties of the first passage time density. The exami- 
nation of the first passage time density shows that both 
the median and the width (defined as the interquantile 
distance) of the first passage time density depend in a 
non-monotonic way on the value of the stability index a. 
There exists such a value of the exponent a characteriz- 
ing the jump length distribution for which the distribu- 
tion width is minimal and also the median is minimal. 
This non-monotonic behavior of the median and the dis- 
tribution width is the consequence of a non-equilibrium 
character of the a-stable Levy type noises. 
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